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We present first data on event-by-event fluctuations in tlie 
average transverse momentum of cfiarged particles produced 
in Pb-I-Pb collisions at the CERN SPS. This measurement 
provides previously unavailable information allowing sensitive 
tests of microscopic and thermodynamic collision models and 
to search for fluctuations expected to occur in the vicinity of 
the predicted QCD phase transition. We find that the ob- 
served variance of the event-by-event average transverse mo- 
mentum is consistent with independent particle production 
modified by the known two-particle correlations due to quan- 
tum statistics and final state interactions and folded with the 
resolution of the NA49 apparatus. For two specific models 
of non-statistical fluctuations in transverse momentum limits 



are derived in terms of fluctuation amplitude. We show that a 
significant part of the parameter space for a model of isospin 
fluctuations predicted as a consequence of chiral symmetry 
restoration in a non-equilibrium scenario is excluded by our 
measurement. 



I. INTRODUCTION 

The ultimate goal in the study of collisions of heavy 
ions at the highest available energies is the production 
and characterization of an extended volume of decon- 
fined quarks and gluons, the quark gluon plasma (QGP) 
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101 . Due to the high muhiphcities in central Pb+Pb col- 
hsions at 158 GeV per Nucleon, recorded in the NA49 
large acceptance spectrometer, a statistically significant 
determination of momentum space distributions and par- 
ticle ratios can be performed for single events, allowing 
for a study of event-by-event fluctuations [g|j3|. In this 
paper we will focus on fluctuations in the average trans- 
verse momentum of individually measured charged par- 
ticles from event to event. 

One expects that the fluctuation patterns are altered in 
the vicinity of the QCD phase transition [^ji) . This con- 
jecture is supported by recent calculations in an effective 
model of the strong interaction [||J^ , which suggest that 
near a tri-critical point in the QCD phase diagram the 
event-by-event fluctuation pattern in transverse momen- 
tum should change significantly. 

A precise measurement of event- by- event fluctuations al- 
lows for a test of the hypothesis of thermal equilibrium 
and the extraction of thermodynamical properties of 
the system in a model comparison. Model studies ||] 
have shown that non-equilibrium models of nuclear colli- 
sions based purely on initial state scattering can be tested 
by measurements of transverse momentum fluctuations. 
Model calculations on transverse momentum fluctuations 
have been performed in many of the commonly used mi- 
croscopic models of nuclear collisions ]9|-[Tl| , in particular 
focussing on the question of how the fluctuations change 
when going from nucleon-nuclcon to nucleus- nucleus col- 
lisions. 

It has also been suggested that for a thermodynamical 
picture of the strongly interacting system formed in the 
collision, the strength of fluctuations is directly related 
to fundamental properties of the system like the specific 
heat (l^Jl^ and matter compressibility A detailed 
discussion of transverse momentum fluctuations in a res- 
onance gas model can be found in |^ . 
One of the most intensely discussed topics related to fluc- 
tuations at the QCD phase transition is the formation of 
so-called disoriented chiral condensates (DCCs) [ p"5|JT6| 
as a consequence of the transient restoration of chiral 
symmetry, which may lead to a production of pions with 
much larger fluctuations of the charged-to-neutral pion 
ratio than expected from Poisson-statistics. The sensitiv- 
ity of our measurement to these fluctuations is discussed. 
NA49 is currently pursuing two different, but comple- 
mentary approaches to the characterization of the single 
events. In the approach presented here we characterize 
the event by global observables like the mean transverse 
momentum of individually detected charged particles in 
the event, averaging over a large interval in momentum 
space. Global quantities in general also include contribu- 
tions from particle correlations occuring at smaller scales, 
i.e. smaller intervals in momentum space. NA49 is also 
studying a system of differential measures of event mor- 
phology which aim at a multiscale characterization of the 
correlation content of single events, which should even- 
tually provide a decomposition of the global fluctuations 
as a function of scale p7| . 



II. EXPERIMENTAL SETUP AND DATA 
SELECTION 

The setup of the NA49 experiment is described in . 
We used a data set of central Pb-|-Pb collisions that where 
selected by a trigger on the energy deposited in the NA49 
forward calorimeter. The trigger accepted only the 5% 
most central events, corresponding to an impact param- 
eter range of 6 < 3.5 fm. The event vertex was recon- 
structed using information from beam position detectors 
and the fit of the measured particle trajectories. Only 
events uniquely reconstructed at the known target posi- 
tion were used. The NA49 large acceptance hadron spec- 
trometer allows the detection of more than 1000 individ- 
ual charged particles for a single central Pb-|-Pb collision. 
In this analysis particles were selected that had a mea- 
sured track length of more than 2 m in one of the two 
Main Time Projection Chambers (MTPC) outside the 
magnectic field and were also observed in at least one 
of the Vertex TPC's inside the superconducting mag- 
nets. We studied particles in a region of 0.005 < pt < 
1.5 GeV/c and rapidity 4 < j/jr < 5.5. A cut on the 
extrapolated impact parameter of the particle track at 
the primary vertex was used to reduce the contribution 
of non- vertex particles originating from weak decays and 
secondary interactions. We estimate that about 60 % of 
such particles are rejected by the vertex cuts. From a full 
simulation of our apparatus using a GEANT ||l^ based 
Monte-Carlo code and a parametrization of the detec- 
tor response we obtained an average reconstruction effi- 
ciency of 90%. The average resolution in transverse mo- 
mentum for the particles used here is around 3 MeV/c, 
dominated by multiple Coulomb scattering. The two- 
track resolution was determined using both the simula- 
tion and a mixed-event technique. For particles selected 
by the track cuts the pair detection efficiency drops from 
around 80% at an average distance in the Main TPC of 
d = 2.5 cm to around 20% at an average distance of 
d ~ 1.5 cm. In table 1 the most important parameters 
of the inclusive and event-by-event distributions of ac- 
cepted particles are summarized. Throughout this paper 
brackets {{x)) will denote averages over events and bars 
(x) will denote inclusive averages over all (accepted) par- 
ticles and all events. 



TABLE I. Measured parameters of the inclusive and 
event-by-event accepted particle distributions. Errors are sta- 
tistical only. 



No. of events 


98426 


(iV) 


270.13 ±0.07 


((A^) - (A>2)i 


23.29 ± 0.05 


Pt 


376.75 ± 0.06 MeV/c 




282.2 ± 0.1 MeV/c 
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III. ANALYSIS AND RESULTS 

For each of the events we characterize the observed par- 
ticle distribution in the acceptance region by calculating 
the mean of the transverse momentum distribution of the 
N accepted particles in the event, 

N 

M{pt)^IIN -Y^PT^. (1) 

i=l 

The resulting distribution of M{pt) is shown in fig. 1. 
The distribution of M{pt) has approximately Gaussian 
shape. No significant excess of 'anomalous' events out- 
side the main distribution is observed. For the variance 
of the M{pt) distribution we get 

y(M(pT))/PT = 4.65 ± 0.01%. 

The biggest contribution to the observed variance is ex- 
pected to come from finite-number statistics. The main 
task in the remainder of the paper will be to extract possi- 
ble non-statistical contributions on top of the trivial sta- 
tistical variation from event to event. A first impression 
of the possible size of non-statistical contributions can be 
obtained by a comparison to the same distribution calcu- 
lated for so-called mixed events (solid line in fig. 1). The 
mixed events were constructed by combining particles 
drawn randomly from different events while reproduc- 
ing the multiplicity distribution of the real events. Only 
one track of any original event was used in a given mixed 
event and no further selection was made regarding the 
impact parameter or multiplicity of the original events. 
By construction the mixed events have the same single- 
particle distributions as the real events, but no internal 
correlations. The variance of the mixed event M[pt) dis- 
tribution is therefore determined by finite number statis- 
tics, giving 

{p^-PT^)'^/{W-VW)) = 4.6%. 

The mixed event distribution resembles, very closely, the 
single event distribution, thus suggesting that large am- 
plitude non-statistical fluctuations are small and/or rare. 
To further quantify and study the deviation of the M{pt) 
distribution a number of methods has been discussed re- 
cently p0[ - p3[[^ In this analysis we follow the approach 
suggested in|0]. We define for every particle i 

Zi = PTi - W- (2) 
For every event we calculate 

N 

Z = Y,z,. (3) 

1=1 

With this definitions we use the following measure to 
quantify the degree of fluctuation in mean transverse mo- 
mentum from event to event: 



One limiting case for this fluctuation measurement is par- 
ticle emission according to a parent distribution that re- 
mains unchanged for all events, i.e. every single event is 
just a random sample of finite multiplicity taken from 
the same parent distribution, ^j,^ was defined such that 
for this case a value of zero is assumed. This value also 
corresponds to the fluctuations for an ideal gas of classi- 
cal particles . 




0.3 0.4 0.5 

M(p.j,) (GeV/c) 

FIG. 1. Event-by-event distribution of the mean transverse 
momentum M{pt) of accepted particles in the event (points). 
For comparison, the solid line shows the M{pt) distribution 
for mixed events. 

Our goal is to detect or exclude fluctuations that are 
compatible with changes, event-to-event, in the parent 
distribution in transverse momentum. Such changes 
would in general lead to values of > 0. ^pj,, as 
we will demonstrate later, is also sensitive to internal 
correlations or anti-correlations of particles within single 
events, which result in > or < 0, respectively. 
In our data set we measure a value of 

= 0.6 ± 1.0 MeV/c, 

compatible with zero. The error was estimated by cal- 
culating <I>py separately on independent subsamples of 
approximately 10000 events each. Before discussing the 
implications of the small value measured for for 
the existence of collective non-statistical fluctuations, we 
first turn to investigating the sensitivity of this measure 
to correlations at small scales, which are known to be 
present in these events. 

Previous studies p5| ] have quantitatively analyzed two- 
particle correlations in relative momentum, with quan- 
tum statistics and Coulomb final state interactions giving 
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the strongest eontributions. We have developed an anal- 
ysis procedure to identify the contribution of these known 
effects, folded with the NA49 experimental response, to 
the observed value of ^pj, . This comparison is done as 
follows. 

1. From the original events we construct a sample 
of mixed events with the same overall multiplic- 
ity distribution as the real events. Without fur- 
ther modifications, the mixed events give a value 



of 



-0.2 ± 0.4 MeV/c, again consistent with 



2. In the second step we model the contributions 
from particle pair correlations at small relative mo- 
menta ('small scales'). The effect of Bose-Einstein 
or Fermi-Dirac statistics on transverse momentum 
fluctuations has been discussed in ||2^ . 

The effects of quantum statistics and final-state 
(Coulomb) interactions partially cancel each other 
and are further diminished by particles emitted 
from long-lived resonances, particles originating 
from weak decays and by including combinations 
of non- identical particles. A detailed evaluation of 
two-particle correlations in NA49 is given in [ p5[ . 
To include the sum of all these effects in the mixed 
events we use a procedure that was described in 
[ pGf . In this procedure the momenta of particles 
are altered pairwise to introduce the desired form of 
the two-particle correlation function, making sure 
that in the mixed events on the average the two- 
particle correlations as a function of relative mo- 
mentum closely match those observed in the data, 
only corrected for the two-track resolution. The 
analysis of the modified mixed events provides an 
estimate of the minimal event-by-event M{pt) fluc- 
tuations that we expect as a consequence of the ob- 
served average two-particle correlations. The con- 
tribution from the two-particle correlation function 
alone is A$py = 5 ± 1.5 MeV/c. 

3. In our data set we observe a slight but statistically 
significant correlation between the multiplicity of 
the event N and the average transverse momentum 
M(j)t)- The correlation is characterised by a linear 
correlation coefhcient of 



{{M{pt)-W)-{N-{N)) 
V{M{pt)) ■ V{N) 



-0.03 ±0.01, 



where V{N) is the variance of the multiplicty dis- 
tribution. We therefore observe a slight decrease 
of M(pt) with increasing multiplicity within the 
central collision data set. Introducing this corre- 
lation in the mixed events gives a negligible con- 
tribution to the width of the M{pt) distribution 
iA<^p^ < 1 MeV/c). 



4. We then apply an experimental filter on each of 
the modified mixed events that simulates the in- 
fluence of the two-track resolution and momentum 
resolution of the NA49 apparatus. While the con- 
tribution from momentum resolution is found to 
be negligible for the range of fluctuations consid- 
ered here, the two-track resolution results in an 
effective anti-correlation between particles in mo- 
mentum space and gives a contribution of A^p^ = 
-4 ±0.5 MeV/c. 

Combining all effects, we find that the observed value 
4>py is compatible with independent particle produc- 
tion: Including both the effects of two-particle corre- 
lations in momentum space and the experimental two- 
track resolution leads to a cancellation resulting in a 
very small net contribution. Any additional contribu- 
tions beyond those mentioned above either have to be 
small or cancel with sufficient accuracy to be compatible 
with = 0.6 MeV/c. 

Here it is worthwhile to note that the effects of two- 
particle momentum correlations and two-track resolu- 
tion, as included in our simulations, both are strongly 
multiplicity-dependent and become negligible for mul- 
tiplicities comparable to those observed in p+p colli- 
sions. This suggests that the physical origin of transverse 
momentum fluctuations as measured using ^p^ changes 
when comparing the value of 0.6 MeV/c for Pb±Pb to 
the preliminary NA49 measurement of 5 ± 1 MeV/c for 
p+p collisions pTj . 

To further study the sensitivity of our measurement 
we have introduced explicit non-statistical fluctuations in 
the mixed event sample and studied the response in ^p^ 
as a function of the parameters controlling the strength 
of these fluctuations. By comparing the value obtained 
for observed in such models with that in the data, 
we can determine the sensitivity of our measurement to 
various kinds of fluctuations and eventually derive limits 
on the amplitude or frequency of occurence for fluctua- 
tions in speciflc models. 

In the first model we examine the sensitivity to non- 
statistical fluctuations introduced by scaling the trans- 
verse momentum for all tracks in a given event by a con- 
stant factor X. The resulting change in the px parent dis- 
tribution from event to event resembles that of an event 
by event change in the inverse slope parameter of an ex- 
ponential transverse momentum distribution. We obtain 
a random number x for each event distributed according 
to 
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exp(-(x~l)V(24„J) (5) 



and multiply the transverse momentum of each particle 
in the event by the same factor x. Here the amplitude 
of the event by event fluctuation is controlled by the pa- 
rameter cT/iuc- Adding these fluctuations to the mixed 
events and neglecting two-particle correlations and two- 
track resolution, increases proportionally to u'ji^^. 
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When adding all effects, we find that for this model a 
fluctuation strength afiuc = 1-2% corresponds to — 
7 MeV/c. Given the definition of ^p^, in Eq. 4 and 
net contributions to $pj, as observed and simulated in 
steps 1 to 4 above we can establish an upper limit on 
cTfiuc < 1-2% at 90% confidence level. 
While a precise limit for specific models can only be set 
using a simulation procedure as outlined above, we typi- 
cally find that for various types of non-statistical fluctua- 
tions our measurement is sensitive when the fluctuations 
lead to an effective non-statistical variation in the mean 
transverse momentum from event to event of about 1%. 

We can also use the simulation procedure to establish 
limits on fluctuations of amplitude (Tfiuc that occur only 
in a fraction F of all events. The resulting exclusion 
plot is shown in fig. 2, where the relative frequency F of 
events exhibiting fluctuations of amplitude afiuc is plot- 
ted versus cr//„c. We see that for F = 1 fluctuations of a 
relative amplitude of cr/;uc > 1.2% are ruled out at 90% 
confidence level, whereas fluctuations occuring in 1% of 
the events can only be ruled out for f7//„c > 10%. 
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FIG. 2. Limit on the amplitude of fluctuations in the pr 
parent distribution as a function of the frequency of events 
showing the fluctuation. 

It is important to note that the measurement of fluc- 
tuations in M{pt) is also relevant for models of processes 
that lead to non-statistical fluctuations localized in trans- 
verse momentum. The most widely discussed example of 
such a process is the formation of disoriented chiral con- 
densates (DCCs) ||l5| , [l6|| , which has been postulated as a 
consequence of the possible restoration of chiral symme- 
try in non-equilibrium scenarios for heavy-ion collisions. 
DCC models predict the formation of domains that even- 
tually emit pions where the ratio / of neutral to all pions 
varies as 



P{f) 



1 



2V7' 



(6) 



The models also suggest that pions emitted from DCC 
domains will be preferentially produced at low transverse 
momenta p^ . This provides for a translation of the 
number-fluctuations predicted by the DCC models into 
Pt fluctuations accessible to our experiment. A limit 
on DCC production has already been set by the WA98 
collaboration |2|], based on a study of fluctuations in 
the relative multiplicties of charged and neutral particles 
near mid-rapidity. 

For comparison purposes we used the same DCC model 
as in |g8[, where the DCC production is characterized by 
the probability F to form a single DCC domain in an 
event and the fraction ^ of pions coming from the DCC. 
We make the additional assumption that the DCC pions 
are produced with pT < p^"'^ — ra.,^. The ratio of neu- 
tral to charged pions was chosen randomly according to 
equ. 1^. The isospin fluctuations of pion production from 
DCCs then lead to multiplicity fluctuations of charged 
pions at low transverse momenta and therefore to non- 
statistical fluctuations in M{pt)- For DCCs occuring in 
every event [F — 1) the fluctuations observed in the data 
rule out DCC sizes of ^ > 3.5 %, which is about a factor 
of 5 smaller than the previous limit set in . This limit 
could be further improved by restricting the analysis to 
the region of small transverse momenta. 



IV. SUMMARY 

In summary, event-by-event fluctuations in the average 
transverse momentum of charged particles in the forward 
hemisphere of central Pb-t-Pb collisions have been mea- 
sured. The distribution of average transverse momentum 
per event M{pt) has an approximately Gaussian shape, 
with no excess of 'anomalous' events falling out of the 
distribution. The fluctuation strength in the data is char- 
acterized by a value of ^p^, = 0.6 ± 1 MeV/c. Using a 
procedure based on mixed events we find that the fiuctu- 
ations in M{pt) from event to event are compatible with 
independent particle production modified by the known 
two-particle correlations due to quantum statistics and 
final state interactions and taking into account the re- 
sponse of the NA49 apparatus, without requiring further 
variations in the transverse momentum parent distribu- 
tion from event to event. 

For a model of non-statistical fiuctuations in average pt 
we use a detailed simulation procedure to determine an 
upper limit on the strength of fluctuations occuring in 
every event of c/iuc < 1.2% at 90% confidence level. 
We also demonstrate that high precision measurements of 
charged particle transverse momentum fluctuations pro- 
vide a sensitive test for models predicting the formation 
of disoriented chiral condensates in heavy-ion collisions. 
Finally, we have provided the first measurement to com- 
pare to predictions of event-by-event transverse momen- 
tum fluctuations in thermodynamical descriptions of the 
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stro igly interacting system produced in Pb+Pb collisions 
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at the SPS. 
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